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ABSTRACT 

We report the detection of radio and X-ray pulsations at a period of 51.6 ms from the X-ray source 
RX/AX J2229. 0+6114 in the error box of the EGRET source 3EG J2227+6122. An ephemeris derived 
from a single ASCA observation and multiple epochs at 1412 MHz from Jodrell Bank indicates steady 
^ ' spin-down with P = 7.83 x 1CU 14 s s -1 . From the measured P and P we derive spin-down power 

E = 2.2 x 10 37 ergs s _1 , magnetic field B p = 2.0 x 10 12 G, and characteristic age P/2P = 10, 460 yr. An 
image from the Chandra X-ray Observatory reveals a point source surrounded by centrally peaked diffuse 
emission that is contained within an incomplete radio shell. We assign the name G106.6+2.9 to this new 
supernova remnant, which is evidently a pulsar wind nebula. For a distance of 3 kpc estimated from 
X-ray absorption, the ratio of X-ray luminosity to spin-down power is « 8 x 10~ 5 , smaller than that of 
most pulsars, but similar to the Vela pulsar. If PSR J2229+6114 is the counterpart of 3EG J2227+6122 
then its efficiency of gamma-ray production, if isotropic, is 0.016 (d/3 kpc) 2 . It obeys an established trend 
of 7-ray efficiency among known 7-ray pulsars which, in combination with the demonstrated absence of 
any other plausible counterpart for 3EG J2227+6122, makes the identification compelling. If confirmed, 
this identification bolsters the pulsar model for unidentified Galactic EGRET sources. 

Subject headings: gamma rays: observations — pulsars: individual (PSR J2229+6114) — supernova 
remnants 



1. INTRODUCTION for RX/AX J2229.0+6114. 



^ ■ 3EG J2227+6122 (Hartman et al. 1999) is one of a 
5h ' number of "unidentified" EGRET sources at low Galactic 



2. CHANDRA IMAGE 



latitude, (£,b) = (106.°6, 2.°9), for which a pulsar origin We observed the X-ray source RX/AX J2229.0+6114 for 

is the hypothesis favored by many authors (e.g., Halpern 17,728 s on 2001 February 14 with the Chandra imaging 

& Ruderman 1993; Kaaret & Cottam 1996; Yadigaroglu CCD array ACIS-I. The target was offset by 0/5 from the 

& Romani 1997). Emission attributed to this source was default ACIS-I pointing position in each direction in order 

also detected by COMPTEL (Iyudin ct al. 1997), and to avoid the inter-CCD gaps. A point source was clearly 

probably by COS B (Wills ct al. 1980). We recently detected, surrounded by diffuse emission with a centrally 

presented an X-ray, radio, and optical study of the error peaked morphology (Figure 1). There is evidently an arc 

circle of 3EG J2227+6122 that revealed only one plausi- of emission that resembles similar structures in the Vela 

ble candidate, a highly polarized, flat-spectrum radio shell nebula (Helfand, Gotthelf, & Halpern 2001), surrounded 

superposed on a compact, non-thermal X-ray source with by fainter diffuse emission that is largely confined within 

power-law photon index T — 1.51 ± 0.14 and with no ob- the boundaries of the radio shell, which has a radius of 

vious optical counterpart (Halpern et al. 2001, hereafter ~ 100". The radial profile of the X-ray flux is shown 

Paper I). We concluded that the most likely interpretation in Figure 2. The central point source accounts for 29% 

was a young pulsar and associated wind nebula (PWN). of the total flux within 100". The total 2-10 keV flux 

The distance was estimated as 3 kpc from the X-ray ab- is 1.3 x 10~ 12 ergs cm -2 s _1 , only slightly smaller than 

sorption column density of (6.3±1.3) x 10 21 cm~ 2 . At that the value of 1.56 x 10~ 12 ergs cm' 2 s _1 that we found in 

distance, the spin-down luminosity required to power the Paper I using ASCA. It is possible that the ASCA GIS, 

X-ray and 7-ray luminosities of 1.7 x 10 33 and 3.7 x 10 35 with its large point-spread function, encompassed some 

ergs s _1 , respectively, would be similar to that of the extended diffuse emission that is not easily recovered in 

younger known EGRET pulsars. In this Letter, we re- this short Chandra observation. In combination with the 

port follow-up observations that confirm a pulsar origin non-thermal spectrum of the X-ray nebula, the morphol- 
ogy indicates a "composite" supernova remnant (to which 
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FIG. 1. — A portion of the Chandra ACIS-I image showing 
RX/AX J2229.0+6114. The image is binned into 1" pixels and 
smoothed with a 3-pixel top-hat filter. It is displayed with square- 
root intensity scaling with the contrast set to highlight faint diffuse 
structures. The linearly-spaced overlaid contours correspond to an 
adaptively smoothed map. Also shown are the locations of X-ray 
point sources detected with > 2<r confidence. 

we assign the name G106.6+2.9), although the radio shell 
is probably a shock between the PWN and the surrounding 
medium rather than a supernova blast wave (see §6). No 
thermal emission has yet been detected from this remnant 
at X-ray or optical wavelengths. 
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FIG. 2. — Radial profile of the Chandra ACIS-I image of 
RX/AX J2229.0+6114. Each radial bin is 8 ACIS pixels wide (3."9). 
The profile is normalized by area and corrected for the instrument 
response across the focal plane, including mirror vignetting and inter- 
CCD chip gaps. The pulsar flux is confined to the central bin, around 
which is clearly evident at least two components of enhanced emis- 
sion, one inside and one outside of 20". The one-sigma error bar is 
shown for each data point. The dot- dashed line represents the esti- 
mated background level. The dashed line denotes the approximate 
outer boundary of the radio shell from Paper I. 




There are eight X-ray sources in the ACIS image that 
can be immediately identified by inspection with bright 
stars on the Palomar Observatory Sky Survey. Two of 
these are present on our optical CCD images of the field 
that were discussed in Paper I. We used them to derive a 
correction of l."5 to the X-ray aspect solution, obtaining 
the precise position for the point source listed in Table 1, 
with an estimated uncertainty radius of 0."5. 

3. RADIO PULSAR 

The position of RX/AX J2229.0+6114 was observed 
with the 76m Lovell radio telescope at Jodrell Bank on 
2001 February 27 and 28 for 1.5 and 2.4 hr, respectively. 
The observations used a cryogenic receiver at a central fre- 
quency of 1412 MHz and a filter bank spectrometer with 
64 channels, each 1MHz wide, for each of two orthogonal 
polarizations. After square-law detection and summing 
of complementary polarizations, the total-power levels in 
each of 64 channels were integrated every 1 ms, 1-bit digi- 
tized, and written to tape for off-line analysis. 

Standard search algorithms detected an unmistakable 
pulsed signal at a barycenter-corrected period P = 
51.6235 ms and dispersion measure DM = 200 ±10 
cm~ 3 pc with signal-to-noise ratio (S/N) of 13.9 and 18.6 
on February 27 and 28, respectively. These S/N values 
are consistent with each other after scaling for the relative 
integration times, suggesting that the pulsar does not scin- 
tillate markedly. Based on this and the known telescope 
sensitivity, we estimate that the period-averaged flux den- 
sity of the pulsar at this frequency is 0.25 mJy, with ~ 30% 
uncertainty. This is a rather low value, fully consistent 
with its non-detection in previous "all-sky" surveys, and 
with the median rms of 0.12mJy reported for the VLA im- 
age in Paper I. The radio pulse profile (Figure 3) consists 
of one featureless peak with full-width at half-maximum 
of 0.08 ± 0.02 cycles. 

We began regular timing observations of this source, and 
a phase-connected solution to timcs-of-arrival spanning 1 1 
days in 2001 March yields the period listed in Table 1, and 
P = (7.80 ± 0.03) x 10~ 14 s s" 1 . The more accurate P 
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FIG. 3. — Radio pulse profile of PSR J2229+6114 at 1412 MHz. 
The instrumental resolution is ta 0.02 of the period. Phase zero is 
arbitrary. 
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listed in Table 1 is derived from a fit to this period in 2001 
March and the period subsequently found in the X-ray 
data obtained 1.6 yr earlier (§4). Thus far, our data do 
not meaningfully constrain the amount of timing noise or 
glitch activity that may be present in PSR J2229+6114. If 
we assume a distance of 3 kpc as estimated from the X-ray 
spectrum in Paper I, the free electron/distance model of 
Taylor & Cordes (1993) would predict DM = 75 cm" 3 pc, 
significantly less than the observed DM = 200. Con- 
versely, for the observed DM the model predicts a dis- 
tance of 12 kpc. It is known that for individual objects 
the model distances/DMs can be wrong by a factor of a 
few, and this may be the case for PSR J2229+6114. In 
fact, the rule of thumb of 0.1 free electrons per hydro- 
gen atom, combined with the measured DM, yields a pre- 
dicted TVh = 6x 10 21 cm~ 2 , in exact agreement with the 
X-ray measured value from the ASCA spectrum. While 
the conversion from column density to distance is itself 
uncertain, we tentatively retain the previous "X-ray" esti- 
mate of d — 3 kpc, while acknowledging that considerable 
uncertainty remains. 
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FIG. 4 — X-ray pulse profile of PSR J2229+6114 in the 0.8-10 
keV band from the ASCA GIS. The instrumental time resolution is 
comparable to the width of one phase bin. Phase zero is arbitrary. 



4. X-RAY PULSAR 

Following the radio pulsar discovery in 2001 February, 
we searched the 113,700 s of ASCA GIS data obtained on 
1999 August 4-7 for pulsations at periods slightly shorter 
than the radio, incorporating in the folding the P that 
would be necessary to connect to the radio period obtained 
1.6 yr later. Using the Z 2 test (Buccheri et al. 1983) to 
gain sensitivity to n harmonics of the fundamental pe- 
riod, we found an unambiguous signal with Z 2 = 99.4 
at P = 51.6196134 ± 0.0000005 ms (MJD 51394.365), re- 
quiring P = 7.827 x 10~ 14 s s _1 to connect to the radio 
ephemcris. The Z 2 statistic is distributed as x 2 with 2n 
degrees of freedom. The probability of a spurious detec- 
tion is 6 x 10 -15 per trial, while only ~ 10 4 independent 
trial periods were searched. Figure 4 shows the folded 
light curve in the 0.8-10 keV band, extracted from a re- 
gion of radius 4' around the source. There are evidently 
two pulses approximately 180° apart in phase, although 
they are not symmetric. The rise time of the higher pulse 
is unresolved; the bin size in Figure 4 (2.55 ms) is compara- 
ble to the instrumental time resolution. (Approximately 
half the data were obtained with 0.5 ms resolution, and 
half with 3.9 ms resolution.) Since most of the extracted 
flux is diffuse, it is not possible to measure the true pulsed 
fraction accurately. After subtracting background from 
an annulus surrounding the nebula, the ASCA GIS light 
curve has a pulsed fraction of 22%, but since the Chandra 
image indicates that only 29% of these photons are from 
the pulsar, the true pulsed fraction must be at least 75%. 

5. OPTICAL OBSERVATIONS 

Optical images of the location of PSR J2229+6114 from 
the MDM 2.4m telescope were shown in Figure 3 of Pa- 
per I. The precise Chandra position of PSR J2229+6114 
is now seen to fall only l."5 north of the R — 17.3 star 
marked A in that figure, and a similar distance south of 
a fainter star of R = 21.3. Given the uncertainty of only 
0."5, the X-ray position is inconsistent with either of these 
stars. Star A is, in fact, a highly reddened A star as we 
determined using the KPNO 2.1m telescope and Goldam 



spectrograph. The proximity to this bright star will make 
any further optical search for the pulsar difficult from the 
ground. Although our images reach a limiting magnitude 
of R — 24.5 in good seeing, a more conservative limit of 
R > 23 probably applies at the location of the pulsar in 
the wings of the bright star. 

6. DISCUSSION 

Applying the standard magnetic dipole model for ro- 
tation powered pulsars, we can use the measured P and 
P of PSR J2229+6114 to derive the spin-down power 
E = Iflfl = 2.2 x 10 37 ergs s _1 , the magnetic field 
B p = 3.2 x 10 19 (PP) 1 / 2 = 2.0 x 10 12 G, and the char- 
acteristic age t = P/2P = 10,460 yr. Compared to 
known 7-ray pulsars PSR J2229+6114 is second only to 
the Crab in spin-down power, and it is significantly more 
luminous than the Vela pulsar (PSR B0833-45), which has 
E = 6.9 x 10 36 erg s" 1 . For a distance of 3 kpc estimated 
from X-ray absorption and radio pulse dispersion, the ra- 
tio of total 2-10 keV X-ray luminosity to spin-down power 
is « 8 x 10~ 5 , smaller than that of most pulsars, but sim- 
ilar to Vela (Hclfand et al. 2001). If PSR J2229+6114 is 
the counterpart of 3EG J2227+6122, then its luminosity 
above 100 MeV is « 3.7 x 10 35 ergs s _1 , and its efficiency 
of gamma-ray production, if isotropic, is 0.016 (<i/3kpc) 2 . 
For Vela, the same analysis gives a 7-ray efficiency of 
0.03 (d/500pc) 2 . 

Among the pulsars that are either reliably or prob- 
ably identified with EGRET sources, there is a trend 
(Thompson et al. 1997, 1999) in which the efficiency of 
> 100 MeV 7-ray production increases with decreasing 
spin-down power (E oc B 2 /P 4 ), or cquivalcntly open field 
line voltage (<& oc B/P 2 ). Their efficiencies range from 
0.002 for the Crab, to of order unity for the middle-aged 
pulsars Geminga and PSR B1055-52. As the source of 
3EG J2227+6122, PSR J2229+6114 would have an effi- 
ciency in accord with the established pattern if its dis- 
tance were close to our estimate of 3 kpc. If, on the other 
hand, the distance were as large as the DM estimate of 
12 kpc, PSR J2229+6114 would have to be more efficient 
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than either the Crab or Vela, but not so much more as 
to rule it out as the source of 3EG J2227+6122. While 
PSR J2229+6114 is apparently similar to the Vela pulsar 
in 7-ray luminosity, its 7-ray spectrum is characterized 
as a power-law of photon index T = 2.24 ± 0.14 (Hart- 
man ct al. 1999). This is steep compared to Vela and 
all other EGRET pulsars except the Crab, for which T = 
2. 19 ±0.02. The X-ray pulse profile of PSR J2229+6114 is 
also reminiscent of that of the Crab, with its two unequal 
peaks. 

Iyudin ct al. (1997) reported the detection of a source 
in the 0.75-3 MeV band with COMPTEL, coincident with 
3EG J2227+6122 but with a much larger error box. This 
detection is consistent in flux with an extrapolation of the 
EGRET spectrum, but it would exceed an extrapolation 
of the 2-10 keV spectrum of RX/AX J2229.0+6114 to 
1 MeV. The soft gamma-ray emission in this region there- 
fore deserves more detailed study. PSR J2229+6114 may 
be one of the brightest pulsars at 1 MeV. 

As a supernova remnant, G106.6+2.9 still has the pecu- 
liar property in the radio of shell morphology, but with an 
extremely flat radio spectrum (a r « 0.0, Paper I). Since 
the X-ray emission appears to be largely confined within 
the radio shell, and since the shell is too small to be the 
blast wave of a 10 4 year old supernova remnant, we con- 
clude as was proposed in Paper I that the radio emission 
comes from a shock driven into the surrounding medium 
cither by the motion of the pulsar or by the expansion of 
the PWN. In the bow-shock interpretation, we can relate 
the now known spin-down power, E = 2.2 x 10 37 ergs s _1 
assumed to be carried almost entirely by the PWN, to the 
velocity of the pulsar v p , the ambient density nn, and the 
radius of the shock r via E = A-k r\ c p . Thus 

E= 2.2xl0 37 (£-) ( * V (_J* crgss- 1 . 

VO.Ol/ V3kpcy VgOkms- 1 / 5 

We measured 1/7 for r , and we assume a relatively low- 
density medium, as might be appropriate at a z-height of 
150 pc or in a cavity previously evacuated by a supernova 
explosion. However, the production of the particle energy 
spectrum needed to explain the flat radio spectrum is still 



not understood theoretically. 

7. CONCLUSIONS AND FURTHER WORK 

We have discovered the likely source of the 7-rays from 
3EG J2227+6122, a young and energetic 51.6 ms pulsar 
that has properties in accord with those of known 7-ray 
pulsars. It may be possible to prove the association by 
finding the corresponding pulsed signal in the EGRET 7- 
ray photons. It is unlikely, however, that the ephemeris of 
such a young pulsar is stable enough to extrapolate back 
with phase coherence to the several epochs of EGRET 
exposure on this source. While we have ongoing radio 
and X-ray observations scheduled to assist in this anal- 
ysis, considerable searching of parameter space will un- 
doubtedly be necessary. Future observations with GLAST 
can definitively test the identification, and can provide 
an excellent pulse profile for comparison with the X-ray 
and radio. Even so, it is important to realize that since, 
as described in Paper I, we have already made a sensi- 
tive search of the entire error circle of 3EG J2227+6122 
without finding any plausible alternative X-ray counter- 
part to much fainter limits, it is more conservative to 
accept the identification with PSR J2229+6114 than to 
doubt it. The only other object that we are aware has 
been considered for identification with 3EG J2227+6122 
is a Be star/X-ray binary with an unknown rotation pe- 
riod (SAX J2239.3+6116 = 4U 2238+60, in 't Zand et al. 
2000), but since that X-ray source is l.°5 from the 7-ray 
centroid, or 3 times the 95% error radius, it is an unlikely 
candidate. Meanwhile, PSR J2229+6114/G106. 6+2.9 is 
sure to become a well-observed example of a young PWN, 
more of which are needed to study the poorly understood 
processes of particle acceleration in pulsars, spin-down 
evolution, MHD flows in PWNe, and interaction with the 
surrounding medium. 

We thank Christine Jordan and Benjamin Collins for 
help with reduction of radio data and X-ray data, respec- 
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ments. This work was supported by NASA grants NAG 5- 
9095 and CXO GO1-2049X to the Columbia Astrophysics 
Laboratory. 
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Table 1 

Parameters of PSR J2229+6114 



Parameter 



Value 



R.A. (J2000) 
Decl. (J2000) 
Galactic Longitude 
Galactic Latitude 
Period (s) 

Period Derivative (s s -1 ) 
Epoch (MJD) 

Dispersion Measure (cm~ 3 pc) 
Distance a (kpc) 

Spin-down Luminosity (ergs s _1 ) 
Characteristic age (yr) 
Magnetic Field (G) 



22 h 29 m 05?28(7) 
+61°14'09."3(5) 
106.°65 
2.°95 

0.05162357393(6) 
7.827(2) x 10~ 14 
51980.0 
200(10) 
- 3 
2.2 x 10 37 
10,460 
2.0 x 10 12 



a Scc text. 



